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^SJ ' A new approach to the quenched propagator in QED beyond the IR hmit 

C ■ 

' is proposed. The method is based on evolution equations in the proper time. 

o ■ 
(N ■ 

: I. INTRODUCTION 

(N . 
> ■ 

in ; 

\0 \ We propose a new non-perturbative analysis of quenched QED beyond the infrared Bloch 

O . 

^ ! Nordsieck approximation using the world line formalism. The main idea is to formulate 

<N ; 

^ , evolution equations in the proper time, inspired by the rather intuitive resulting classical 
, picture. As a warm up we obtain in a simple way the renormalized quenched propagator 
D \ in the infrared (IR) limit. From there we proceed to improve on this IR approximation, 
^ I which we achieve in the form of integro-differential equations, our main results. They are 
^ I of course difficult to solve, but they have an appealing physical interpretation and suggest 
a simplified numerical analysis. 

The paper is organized as follows. In section || we study the infrared behavior of lepton 
propagators for several cases. First we consider the Bloch Nordsieck approximation as given 
by uniform proper time paths, which exhibit the classical image in the worldline. This 
suggests the use of the proper time evolution equations as a short cut to go beyond the 
infrared approximation. The idea is carried out for scalar QED and then for the full Dirac 
algebra. The lesson of this section is a better understanding of renormalization and the 
fact that the scalar case contains the basic physics, which will be used in the next section. 



section III where, we develop further the proper time formalism to go beyond the infrared 



1 



limit. The results are summarized in section 



devoted to the conclusions. 



II. INFRARED BEHAVIOR: THE CLASSICAL IMAGE OF THE WORLDLINE 

FORMALISM 

In this section we will study the classical Block-Nordsieck approximation and their ex- 
tensions for scalar and fermionic QED. Our main point is given in Equation (|12|), where an 
evolution equation is introduced for the quenched propagator. Its solution is the well-known 
Block-Nordsieck result. The extensions to the scalar and full Dirac algebra are carried out 
without difficulty. 



A. Bloch Nordsieck approximation 

The Bloch-Nordsieck model has been repeatedly used to study the infrared behavior of 
QED. Here we will consider it from the path integral point of view as a convenient starting 
point. The model has been extensively studied in the literature 

In the Bloch-Nordsieck model one substitutes the Dirac's gamma matrices by a constant 
vector 7^ —>■ u^. This constant vector will be identified to the (constant) velocity of the 
electron, see below. With this substitution the model is completely integrable and the elec- 
tron Green function can be cast in terms of elementary functions. The electron propagator 
in the Bloch-Nordsieck approximation is given by 

d 

[u^,{i-Q^ + eAf,{x)) -m]G{x,y) = -5{x-y), (1) 

where v4^(x) is the electromagnetic potential. We will solve this model by introducing a 
proper time (Schwinger's method [|]) 0J^,|,3. Define 

POO 

G{x,y) = i dTG{T,x,y), (2) 
Jo 

then G{t, x, y) verifies the equation 
dG d 

~ = + eA^W) - m]G(r, x, y) = HG{t, x, y), (3) 



with initial conditions 

G(T = 0,x,y) = 5(x-y). (4) 

The formal solution is given by 

G{T,x,y)^ex.p{-iTn} S{x-y). (5) 

Here 7i is clearly seen as a Hamiltonian which gives the evolution in the proper time of the 
Schrodinger-like equation. This formal expression can be cast into a path integral in the 
proper time. We can write 

G{T,x,y) = N J Dx{t)Dp{t) cxp{i J dT{px — T-C{x,p))} 

— ^ j Dx{t)Dp{t) exp{i j dT{px — up — euA(x) + m)}. (6) 

The electromagnetic field is coupled to the current j^(^) — Jq drS^z — x{t)) and can be 
integrated out, since the resulting integral is gaussian 

< G{t, X, y) >= j DA^ G(t, x, y) exp{-z j dx-^F^ + ^{dAf}, (7) 

where A is a gauge-fixing parameter and the brackets <> mean the average over the gauge 
fields. 

The integration of the A fields gives a non-local lagrangian for < G > 



G{T,x,y) > =N j DxDpexp{i j dT{px — pu + m)} 



Dni,{k) is the photon propagator in the gauge A. 

The integral over p gives a delta function, 5{x^{t) —u^)), which implies that the particle 
moves with constant velocity, u^, as stated previously. We see that the electron in the Bloch- 
Nordsieck approximation behaves as a classical particle, the only path that contributes to 
the propagator is the classical path and quantum fluctuations are exactly canceled. This 
fact is most easily seen in the above formalism. The triviality of the quantum corrections 



is due to the appearance of the delta function, which in turn is related to the triviality of 
the classical hamiltonian. For the BN model, the classical path and the uniform velocity 
path coincide, independently of the external potential, Afj_{x). Alternatively one may notice 
that the Eq.(|^) is first order on the variable x. As is well known from the theory of partial 
differential equations, the solution admits a particle interpretation. This is not true for 
the cases considered below. The uniform path is only an approximation valid when the 
momentum interchanged with the gauge fields is sufficiently small such that the path can 
be considered uniform, i.e. in the infrared domain. 

Therefore the non-local term must be evaluated in the classical path alone. It gives 

N.L. = r dndT2-. (9) 

After renormalization we find 

< G{t,p) >= ir"" exp{—iT{m — up)}, (10) 

where a = (3 — A)/(27r)Q:. Integrating over r, we obtain finally 

poo i 

< Gip) >= ^l dr< Gir,p) >= r(l + a) ^^_^^^^^^ . (11) 

This is the well-known Bloch-Nordsieck result. In the Yennie-Suura gauge, A = 3, the elec- 
tron propagator reduces to the free case. Notice that < G{p, r) > verifies the homogeneous 
evolution equation 

- i^^^klll^ = (up -m)< Gip, r) > +- < Gip, r) >, (12) 
or T 

which has of course the Eq. (p!oD above as a solution. This fact will be used later on to go 
beyond the infrared approximation. 



B. Scalar QED 



Before considering the spinor electrodynamics, let's study the somewhat simpler case of 
scalar QED. We will study in detail the renormalization for this case. 

4 



The scalar propagator, coupled to the electromagnetic current, is given by 

{{d, - eA^{x)y - m^)G{x, y) = 5{x - y), (13) 
as before we introduce the Schwinger proper time 

.dG{^^ = {{d, - eA,{x)f - m')G{x, y), (14) 
G{t, X, y) has the solution 

G(t, x, y) = exp{— iTir} (5(a; — y) (15) 

which can be expressed as a path integral 

G(T,x,y) J Dx{r)Dp{r)exp{i J dr{px + n{x,p))} 

^ N J Dx{r)Dp{T)exp{i J dT{px + {p + eA{x)f - 171"^)}. (16) 

Now p appears quadratically and the integral over it is no longer a delta function. As 
expected, the classical path is not the only contribution to the path integral and quantum 
corrections are relevant. However, for the moment, we will consider the contribution given 
by the uniform path alone, as the main contribution in the infrared limit. We first shift the 
momentum integral 11 = p + eA[x) 

G{r, x,y) = N J Dx{t)DU{t) exp{i £ dT{Ux + -m^ - exA^{x))}. (17) 

and as before, the electromagnetic field is coupled to the current jn{z) = J dTX^{T)5{x{T) — 
z). The average over gauge fields gives therefore, the same non local term 



N.L. = / drdr' j j^x,{t)x.{t')^^ e^'^<^^-<^'^\ (18) 



Now, we make the uniform path approximation, x^(t) — u^ — const, which gives 

dk 



N.L. 

Integrating over r and r' gives 



The integral is IR finite, at low momentum the divergent denominator cancels out with the 
1 — cos(/cnr) in the numerator. We will renormalize the ultraviolet divergent k integral using 
dimensional regularization, expanding around four dimensions, D = 4 — e, we obtain |^ 

N.L. = (21) 

where 

^ = - - ^(3 - fl)) = A„ + + . . . . (22) 

and Aq = — l/(47r^). The scalar propagator can be written 

< G{x, y) >= [ -J^e'P^^'-y^ r dr exp -iT{m' - exp -e^-{^\y. (23) 
J [271)^ Jo e 

Expanding the exponential in powers of e we obtain 

n=0 

Each term in the series can be integrated and gives 



< G{p, t) >= / dr exp -iT{m^ - p") V (^_)"(;iV)"^ (24) 
-'0 rr^, e 



< G{p^^) >= E + ne)(4^)"^ (25) 

(m^ — p^) n!e" — 

The n-th term can be written as the n-th power of 

Expanding now in powers of e we arrive at 

-e'(^^ + A,- Aoj + Ao log((2z^2)/(^2 _ ^2)) ^ ^^^^^^ ^ (27) 

and we can choose the counterterms in such a way to cancel the term Ao/e + Ai — Aqj. The 
scalar propagator after renormalization is given by 

< Gnip) >= ' e-aiog(MV(r»^-P^))^ (28) 
— 

which coincides with the Block-Nordsieck result, as expected. Notice that care has to be 
taken in the order of the limit e — > and the integral over r. We will see next that including 
the Dirac algebra is under control, this will allow us to go back to the scalar case for the 
purposes of the present article. 



C. Full Dirac algebra 

Including the full Dirac algebra is far from trivial. It has been done in references 

We will follow the method of Ref. [|l^. The Dirac equation is 

- e-i^A^{x) - m)G{x, y) = 6{x - y). (29) 

Introducing the proper time as before, (see Eq.(6) from ref. |1T0[| ) 

G{x,y)= f/Ty dxexp{--(m2r + mx)}exp{-((7n)2r + 7nx)}, (30) 

where II'^ = — eA^{x) and x is a Grassman variable. Introducing a set of auxiliary 
Grassman variables we can write a path integral for the propagator 

G(x,i/) = exp{7^} dT JdxJ DxDpD^exp{z dripx + - H) - m^T)}, (31) 

where TC = ^(—11^ + 2ieFfj,y^^^'' — |;n^^^x). The verifies the boundary conditions 

■6^(0) +^/i(^) = r^. Shifting the momentum as before and integrating the resulting gaussian 

integral over 11 we obtain 

d fir 
Gix,y)=exp{^—} dT J dxJ^iT) exp{--im^T + mx)} J DxD^ exp{tS[x, ^]}, (32) 

where ^] is given by 

S[x,^] = J^dr{-^x' + - ^eF,,eC - ^i^x) - ^e(0)e(^). (33) 

As before the gauge fields appear linearly only and the integration over them can be done 
exactly. 

< G{x,y) >= exp{7^} j^'rfT J dxM{T) J DxD^ exp{zS,s}, (34) 
where the action Scs is given by 



SeS= I dT{--{xf + t^i) 

Jo 2 



(2 



— x{ti) and so on. 

We will calculate the action for the uniform path, as done before. 

From the classical equations of movement for the free case (e = 0) we get that the 
uniform path in our case is given by 

Ur) = U^) - ^P,X- (36) 
Introducing this solution in the path integral and after some algebra we arrive at 

< G(x,y) > =exp{7^} J ^^e'Pi^-v) J^"^ dT J dxexp{^{m'T + mx)} 

eM^P'T + ipmx - e'^i^iTpYil + ^^e)}, (37) 
where A is given as before Integrating over x ^-nd applying the d/dV we arrive at 



< 



dp An(x-v) r .rr. (-\...2rr.^^ r ^.2^ .2^, 

(27r)^ 

(Y(^ + P7) + ^(/^We^). (38) 

If we set e = we obviously reproduce the free fermion propagator. For e 7^ care must be 
taken in the regularization procedure since the integration over T and the limit e — > do 
not commute. If we boldly take the limit first and integrate we arrive at 

<''W>-Z(^H^r(i + „)-^(™ + p,(i + ?^^)), (39) 

where Z is the renormalization constant Z = exp(— e^/e74 + . . .) 
Taking correctly the limit after the integration we arrive at 

< G{p) >= Z i^^^" , r(l + a)^^ (m +P7) , (40) 

which is again the Block-Nordsieck result and reproduces the IR logs upon expansion on 
powers of a. 
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III. BEYOND INFRARED BEHAVIOR: PROPER TIME EVOLUTION 

EQUATIONS 



We will now write down the evolution equations for a scalar coupled to a scalar field. We 
will use this simpler case to demonstrate that the evolution equations allow us to go beyond 
the infrared domain. The obvious next step would be to construct evolution equations whose 
perturbative expansion give exact results up to a given order. Let's start by writing down 
equations valid up to first order. Further generalizations to the spinor case and higher orders 
will be left for future work. 

As before we start with the scalar propagator as given in Eq.(^) in the proper time 
formalism, but now for a scalar field coupled to a scalar field Ti = — iv? + g(f){x). The 
propagator verifies the evolution equation 

-'^ = [(^2 - + 9<P(^)P{r, X, y), (41) 

where g is the coupling constant to the scalar field We want to average over the field 

(p with a weight given by the free action. 

<0>= J d(l)Oexp{i J dxC^}, (42) 

= 0(9^ — 'rn?^)(j). Therefore we get 
.d<G> 



dr 



{dl -m')<G>+g< 0(x)G(r, x,y)>. (43) 



In the infrared limit the average < (j){x)G{T,x,y) > is simply given by a/rG, as shown 
above by Eg . (p^ . But in general it can obviously not be written as a simple operator acting 
on G. So one has to estimate the value of < (f){x)G{T,x,y) > by introducing as before a 
path integral representation for G. 

G{T,x,y) = N J DxDpexp{—i J drpx + p"^ — + g(p{x)}, (44) 

then 
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< 



(f){x)G{T,x,y) > — J D(f)DxDp(f)(x) exp{—i J drpx + — + g(l){x)} 

exp{i j dx(f){x){p'^ — m^)(f){x)}, (45) 

introducing the current jo{z) = g J drd^z — x{t)) the average can be written 

< (j)(x)G(T,x,y) > = [ DxDpexp{—i [ drpx+p^ — m^} 

oj{x) J Jo 

J D(f) exp{i J dx(f>{p^ - ml)(j) + (46) 

where the derivative is evaluated at j = jo- The integral over 4> is gaussian an after per- 
forming the functional derivative we arrive at 

< (f){x)G{T,x,y) > — j DxDpexp{—i j drpx + p^ — m^} 

exp{— j dTidT3D{x{Ti),x{T2))}g dTQD{x,x{TQ)), (47) 

where D{x, y) is the free scalar propagator of field 0. This expression can be simplified fur- 
ther if we make some approximations and constitutes the starting point to obtain evolution 
equations. To start with one can neglect the non-local term in the exponential since it is 
order g'^. Then it is straightforward to evaluate the functional integral since all integrals are 
gaussian. After some algebra we arrive at 
< 0(x)G(t, x,y) >^ g j dro J dkoD{ko)e'''°'' J DxDpexp{-i J drp^ + px - + jx}, 

(48) 

where j(r) = —koS{T — tq). The integral over x can be done and gives now a delta function 
S{p—j), i.e. the path is still uniform but now presents a jump at r = Tq where the momentum 
changes from p to p — ko. The electron emits a hard photon of momentum ko. Evaluating 
the action at this specific path gives the total contribution, in this approximation 

< (f>{x)G{T, x,y)> ^g£ dro J dkoDiko)e'^''°-P^^^-y^ 

exp{-i((p - A;o)^(t - Tq) - m'^{T - Tq) +pVo - mVo)}. (49) 

We have separated the contributions from r — Tq and Tq to make it clear that they are both 
the free scalar propagators , i.e. 
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<<P{x)G{T,x,y) >=gJ^dro J dA;oL'(fco)e^'°(^-^)Go(ro,p)Go(r - ro,p - M- (50) 
This result is valid up to first order on g. Then the evolution equation, vahd up to this order 



IS 



.d<G> 
i- 



(p^ - m^) < G > +^2^ dTQ J dkoD{ko)Go{To,p)GoiT -To,p- ko). (51) 



dr 

The solution of this equation up to order g'^ gives the exact one loop scalar propagator, as 
it should. A natural extension of this result would be to change the free propagators. Go, 
in the integral by the exact ones G with the hope that the main higher order contributions 
will be included in this way. We arrive, therefore, to our final evolution equation 

- i — — = {p^ -m^) <G> +g^ dr^ j dkoD{ko)G{ro,p)G{T -To,p- ko). (52) 

The solution of this equation gives the exact result up to first order. We think that higher 
orders are approximately taken into account. In the infrared limit, the equation reduces to 
the Block-Nordsieck case. 



IV. CONCLUSIONS 

We have arrived at an expression for the quenched propagator beyond the infrared do- 
main. It has a clear physical interpretation as one can see by looking at the diagrams 
resummed by the above equation. From direct inspection, one sees that the vertex function 
is evaluated exactly at the one loop level, and that there are an infinity number of such hard 
loops. Therefore, hard photons do not connect propagators with different hard vertex, i.e. 
we have some kind of rainbow expansion, where overlapping loops are not included, as shown 
by Fig.|l[ But there are in addition an infinity number of infrared photons which join any 
two points in the propagator and make the infrared limit exact. Since the proper time in the 
diagram is related to the virtuality of the propagator and we have, in the above equation, 
an integral from to the upper limit r, we expect that higher order loops are calculated 
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only up to some minimum virtuality, i.e. one expects an ordering in the virtuality of the 
diagrams that contribute to the above equation. 





FIG. 1. Diagrams included (upper graph) and not included (lower graph) on the evolution 
equation 

New insights in the nonperturbative regime are always welcome, even if quenched. In 
fact, better understanding of quenched approximations are essential for lattice simulations. 
The non-positron approximation is also useful in analytic attempts to fundamental problems, 
like stability of relativistic QED, where renormalization is one of the main problems [jl2 



Our result in Eq.(0) can have many practical applications like a systematic way to calculate 
hard loop corrections with exponentiation of soft photons. They will depend on our ability 
to solve or approximate analytically the above equation. Of course Eq.(|52D constitutes a 
new approach and can serve in any case as a starting point for new numerical analysis. 
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